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Abstract 

We present a formulation of domain-wall fermions in the Schrodinger func- 
tional by following a universality argument. To examine the formulation, we 
numerically investigate the spectrum of the free operator and perform a one- 
loop analysis to confirm universality and renormalizability. We also study the 
breaking of the Ginsparg- Wilson relation to understand the structure of chiral 
symmetry breaking from two sources: The bulk and boundary. Furthermore, we 
discuss the lattice artifacts of the step scaling function by comparing with other 
fermion discretizations. 

1 Introduction 

In the study of CP violation by CKM unitary triangle analysis, hadron ma- 
trix elements of four- fermion operators, such as i?Ki play a vital role. Accurate 
calculations of this quantity from first principles are an important task for the 
lattice QCD community. In such calculations, having chiral symmetry is cru- 
cial to avoid an operator mixing problem which causes uncontrollable systematic 
errors. Although lattice chiral fermions [H [21 13] are a clean formulation, they 
require enormous computing power to perform dynamical simulations. In com- 
parison, ordinary fermion formulations, like Wilson type fermions and starggared 
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fermions are relatively cheap. Nowadays, however, thanks to the development 
of computer architecture and algorithms, dynamical simulations with lattice chi- 
ral fermions have become feasible even for three flavors [4J. In particular, the 
RBC/UKQCD collaboration |5] is currently using domain- wall fermions (DWFs) 
to compute Sk- In the course of their computation, there are many sources of 
systematic errors which one has to control. Among them, the non-perturbative 
renormalization (NPR) could be serious. At the moment, the collaboration has 
been using conventional schemes, such as, the RI/MOM scheme and its variants 
[HI [7]. However, these schemes potentially contain "large scale problem" which 
requires a quite large lattice volume. To avoid such difficulties, a new scheme 
was invented, known as the Schrodinger functional (SF) scheme [8]. This scheme 
provides a reliable way of estimating errors in the NPR. If one wants to use this 
scheme for the renormalization of Bk given by the RBC collaboration, first of 
all, one has to formulate DWF in the SF setup. This is the purpose of this paper. 

While chiral fermions are useful for computing the bare i?K to avoid the mix- 
ing problem, a formulation for such fermions in the SF setup was a non-trivial 
task because SF boundary conditions break chiral symmetry explicitly. We will 
address this issue in the next section. However, Taniguchi ^ made the first 
attempt to formulate overlap fermions by using an orbifolding technique. Sub- 
sequently he provided a formulation for domain-wall fermions and then he and 
his collaborators ^Oj calculated a renormalized Bk in quenched QCD. Sint [11] 
developed such techniques by combining with a flavor twisting trick. However, 
these orbifolding formulations are constrained by the requirement that the num- 
ber of flavors be even. Thus, apparently such formulations are incompatible 
with current trends toward dynamical three flavor simulations. To overcome this 
difficulty, Liischer [12] gave a completely different approach relying on a uni- 
versality argument, dimensional power counting and symmetry considerations. 
Some perturbative calculations were performed in Ref. [1^. A crucial property 
of this formulation is that there is no restriction on the number of flavors. Since 
only overlap fermions were considered in Ref. |12j . our main purpose here is to 
formulate the other chiral fermions, namely, domain-wall fermions. 

The rest of the paper is organized as follows. Section [2] gives the formulation 
of domain- wall fermions in the SF setup, after a brief review of the universality 
argument. We present several pieces of numerical evidence in Section [3] and 
mto show that our formulation is working properly. We also discuss the lattice 
artifacts for the step scaling function in Section [5l In the last section, we conclude 
by giving some remarks and outlook. 
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2 Formulation 



In the following, we assume that the reader is familiar with the SF in QCD [8|ll4j. 
After giving a brief reminder of the universality argument, we give a formulation 
for DWF and finally check the chiral symmetry breaking structure numerically. 

2.1 Universality argument 

In the massless continuum theory, the Dirac operator D satisfies the anti-commutation 
relation with 75 

-f^D + D-f5 = 0. (1) 
The above is true even in the SF setup, although the boundary conditions, 

P+ij{x) = at xo = 0, (2) 
P_V'(x) = at xo = T, (3) 

with P-i- = (1 ±70)72, break chiral symmetry explicitly. Eq.([T]) means that the 
operator itself does not know about boundary conditions. In the continuum the- 
ory, information such as boundary conditions is embedded in the Hilbert space. 
In fact, the corresponding propagator, which is a solution of the inhomogeneous 
equation, 

DS{x,y) = 6{x-y), (4) 
fails to satisfy the anti-commutation relation. Instead, it follows 

75S'(x,y) + S'(x,y)75 = 

/ d^zS{x,zh^P-S{z,y)+ [ d^zS{x,z)j5P+S{z,y). (5) 

Jzo=0 Jzo=T 

This can be derived by using partial integration on the SF manifold which has 
two boundaries at time slice xq = and T. The non-vanishing right-hand side in 
eq.© shows an explicit chiral symmetry breaking. Since such a breaking term 
is supported only on the time boundaries, the chiral symmetry is preserved in a 
bulk. 

If someone naively tries to formulate chiral fermions on the lattice, one may 
define an overlap operator, for example, with the Wilson kernel in the SF setup 
|14j . However such an operator immediately satisfies the Ginsberg- Wilson rela- 
tion and thus cannot reproduce eq.(l5|) in the continuum limit. This indicates 
that such naive formulation does not work and furthermore may belong to an- 
other boundary universality class which is not what we want. In this way, it is a 
non-trivial task to formulate chiral fermions in the SF setup. 

Some years ago, Liischer proposed a clever way to overcome this situ- 
ation. First, consider the relation for the propagator in eq.©. This indicates 
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that the GW relation has to be modified by boundary effects. Thus one has to 
find a modified overlap operator which breaks the GW relation near the time 
boundaries and correctly reproduces eq.dS]) in the continuum limit. Actually, 
finding such a modified operator is not so hard. However, a new question natu- 
rally arising is how the SF boundary conditions emerge. For the Wilson fermion 
case [H], because there is a transfer matrix, it is natural for fermion fields to 
follow the SF boundary conditions. However for chiral fermions, there is no such 
transfer matrix which can be defined from nearest neighbor interaction in the 
time direction. Therefore it is not an easy task. 

Liischer [12] gave another point of view to see how fields respect the boundary 
condition. In the quantum field theory, the correlation function can tell you what 
kinds of boundary conditions are imposed. As an example, let us see how the 
boundary conditions emerge for Wilson fermions whose action is given by 



x)D^{m)'4}{x) 



(6) 



(7) 



where and V* are forward and backward covariant difference operators re- 
spectively, 

1 



[U{x, p)^{x + a[i) - i){x)] , 
\i^{x) — U{x — afi, fj,)~^ip{x — a/i)] 



(8) 
(9) 



In the SF setup, the sum over x in the action is a little bit subtle. We assume 
that the dynamical fields are 'ip{x) with a < xq < T — a and the fields ^/'(x) with 
xo < and T < xq are set to zero. For this setup, the propagator may be defined 

by 

-4 s: 
a Or 



{r]{x)i){y)) 
r?(x) 



5'4){x) 



For 2a < Xq <T — 2a, eq. ffTTl) turns out to be 

ri(x) = D„{m)ilj{x) 
On the other hand, at xq = a, we obtain 



(10) 
(11) 

(12) 



r/(x) 



1 



-P+tpix) - VoP-il^ix) 



+ 



ip{x) + m'ip{x). 



(13) 
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By substituting eq.([T3]) into eq.p^ with x ^ y we obtain 

ip+(V'(x)v;(y))Uo=a - VoP-{^{x)ij{y))\„ + ■■. = 0. (14) 
In the continuum hmit, the first term is dominant 

ip+(V'(x)V;(y))U=o = 0. (15) 

This shows that in the naive continuum hmit, the Dirichlet type boundary con- 
dition (P+'i/'Uo=o = 0) is stable against the Neumann one (VoP-'(/'|a;o=o = 0), 
and in the end the SF boundary conditions in eq.([2]) emerge. It is plausible that 
similar things happen also for the chiral fermions case, as long as the locality 
and symmetry are kept in a proper way, although we expect that the coefficient 
of the lowest dimensional operators (^P+i/') may be different from the above 
case, and more higher dimensional terms may appear in eq. ()13p . The important 
point here is that continuum SF boundary conditions emerge dynamically in the 
continuum limit of the correlation function. This boundary condition is natural 
and automatically guaranteed to emerge from the dimensional order counting 
argument. Therefore, when we construct chiral fermions in the SF, we only have 
to prepare a modified operator by introducing an additional term which breaks 
the chiral symmetry near the time boundaries. Once this is fulfilled, then such 
an operator automatically turns out to be the desired one in the continuum limit 
without fine tuning. A final important note is that the form of the boundary 
term is irrelevant as long as it will go into a preferred boundary universality 
class. Therefore, there is a large amount of freedom when choosing boundary 
terms and one can use this freedom for practical purposes. 

Following these guiding principles, Liischer [T2] proposed the operator: 

aD^ = 1-]^{U + U), (16) 

U = A(AU + caP)"^/^ [7 = 75C/t^5, (17) 
A = l + s-aP)„(0), a = a/(l + s), (18) 

with the parameter in the range \s\ < 1. D^{0) is the massless Wilson operator 
in the SF. The key point here is the presence of the P term in the inverse square 
root which is given by 

aP{x,y) = 5^^y5xQ,yo{5xo,aP- + 5xQ,T-aP+)- (19) 

Note that this term is supported near the time boundaries and thus called a 
boundary operator. The presence of this term breaks the GW relation explicitly 
and the breaking is given by 

Ab = 75-Dn + -Dn75 - aP>N75-DN- (20) 
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It was shown in ref. that this term is local and supported in the vicinity of 
the boundaries up to the exponentially small tails. 

Although this operator breaks chiral symmetry explicitly, other symmetries 
(the discrete rotational symmetries, C, P and T, flavor symmetry and so on) 
have to be maintained since the boundary conditions in eq. (|2|3p are invariant 
under these symmetries. In addition, this operator has 75-Hermiticity. In this 
way, the universality formulation can avoid breaking important symmetries, such 
as the flavor symmetry. This is a distinctive feature of this formulation compared 
with the orbifolding technique, where flavor symmetries cannot be maintained 
or, there is a constraint on the number of flavors. 

Before leaving this subsection, let us summarized the guiding principles of 
formulating chiral fermions in the SF setup. What we learned from this con- 
struction is that, for an original chiral fermion operator, one has to introduce an 
additional term to break the chiral symmetry and then demand that such break- 
ing only appears near the time boundaries. Furthermore, one must maintain 
important symmetries as well as 75-IIermiticity. Once these conditions are ful- 
filled, it is automatically guaranteed that the such a lattice operator will correctly 
reproduce the continuum results according to the universality argument. 



2.2 Formulation of domain- wall fermions 

Let us apply the guiding principles given in the previous subsection to domain- 
wall fermions. We propose a massles^ domain- wall fermion action 



S = a^^ ^ 'iI){x,s){Dd^f)^s,x's''^{x' ,s'), (21) 

x,x' s,s'=l 

where a massless operator with = 6 for examplcH in four dimensional block 
form is given by 



with the chiral projections. 
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-Pl + cB 














-Pr - cB 


aD^ 


-Pl 








-cB 





-Pr 


aD^ 


-Pl 


-cB 











-Pr 





(22) 



P, 



R/L 



(l±75)/2. 



(23) 



^The mass term can be introduced in the usual way, namely a'*mf ^al^i^' ^)Pr^{^^ ^s) + 

-0(x, Ls)PLlpix, 1)]. 

^We restrict ourselves to an even number of L^, which is the case usually implemented. 
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We also assume that the dynamical fields are ip{x,s) with a < xq <T — a. The 
block elements in eq. (i22|l are four dimensional operators and aD^ is given by 

aZ)w = aD^i-m^) + 1. (24) 

The domain- wall height parameter usually takes a value in a range < am^ < 2. 

An important ingredient here is the presence of B in eq. (j22p terms in the 
cross diagonal elements. The reason for this s-dependence is to break the chiral 
symmetry in a similar way to the usual mass term [15]. As mentioned before, 
such chiral symmetry breaking should be present only near time boundaries, 
therefore, we chose the B term as 

B{x,y) = (^x,y<5zo,j/o75(5xo,a^'- + 5xo,T-a^ + ), (25) 

which is supported near the boundaries. In this way, the time dependence is 
fixed. The spinor structure (75i-*±) is determined by imposing the discrete sym- 
metries C, P and T and Fs-Hermiticity. These requirements are not so strong to 
determine the spinor structure completely and therefore there is some freedom. 
The structure proposed here is only one of many solutions. Actually, we exam- 
ined several choices of the spinor structure in the boundary term and confirmed 
numerically the universal results in the continuum limit for the lowest eigenvalue. 
In the following, we take this boundary term in eq. (|25p . 

The boundary coefficient c is supposed to be non-zero to correctly reproduce 
the continuum theory as we will see in Section [3l It also plays an important role 
to cancel boundary 0(a) cutoff effects and has a perturbative expansion 

c = c(°)+c(i)5o' + 0(<7o')- (26) 

In the same way as in |13j . we tune the first coefficient c^^^ as a function of the 
domain- wall hight ams, 

cW = 0.5089 - 0.0067(am5 - 1) + 0.0488(am5 - 1)^ 

-0.0216(am5 - 1)^ + 0.0673(am5 - 1)"^. (27) 

This is valid in the region where bulk 0(a) can be neg lectedl ,that is, for sufficient 
large L^. In the process of this determination, one needs to define the operators 
of the axial vector current and pseudo scalar density. We give their definition 
together with that of the conserved axial current in appendix IA.3I 



7 



Ls=4 L3=12 1^=24 




Figure 1: The color ranging [—14,0] corresponds to the value of In ||A'^-^''^(a;o,?/o)||spm 
for the zero spatial momentum configuration with the parameters T/a = 20, mi = 0, 
ams = 1 and c = 1. 



2.3 Structure of chiral symmetry breaking at tree 
level 

In this subsection, let us check an important properties of the operator defined 
in the previous section. A reader may worry that even though the additional 
boundary term is localized to time boundary, after integrating over the fifth di- 
mensional degree of freedom such breaking effects may leak into the 4-dim bulk 
and ruin the bulk chiral symmetry. To settle this question, we numerically inves- 
tigate the structure of the chiral symmetry breaking by looking at the breaking 
of the GW relation 

A^"-^^ = 75i^S^^ + i?S^^5 - 2aDi^^\,Dl^'\ (28) 
with the effective four dimensional operator [T6l \T7\ 

det d[^'^ = det[L>DWF/i^Pv]. (29) 

The Pauli-Villars (PV) operator is defined as the massive DWF operator with 
anif = — 1. To obtain the effective operator, first of all, we have to define physical 

^ Actually, we observe that 0(a) improvement program does not work for small Lg and values of 
arris which are far from 1. For example 0(a) terms in Ja ad fp (defined in appendix lA.Sp at tree 
level do not vanish simultaneously with the same value of c^"-* . 
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quark fields 

q{x) = PLij{x,l)+PRi^{x,Ls), (30) 
q{x) = i^{x,l)PR + ^{x,Ls)PL. (31) 

In terms of the propagator of domain-wall fermions defined from 

DnwFSuwFix, y, s, t) = a~^6^^y5s,t, (32) 

that of the physical field is expressed 

[Q{x)q{y)]F = 5'q(x,y) 

= -Pl'S'dwf(2;, y; 1, Ls)Pl + PlSdwf{x, y; 1, 1)Pr 

+ PRSBWFix,y;Ls,Ls)PL + PRSDWFix,y;Ls,l)PR. (33) 

In terms of Sq the effective operator is given by 

aL»S'^ = (l + «'5q)"'- (34) 

In the SF setup, A^^'-* in eq. ()28p contains not only the bulk chiral symmetry 
breaking but also the boundary breaking. The former is supposed to be removed 
by taking Lg to infinity. In such limit, boundary breaking effects remain and 
they are expected to be localized near time boundaries. To see this situation, we 
numerically compute A^^") for a free operator. In the free case, we can perform 
the Fourier transformation for spatial directions. We study the momentum con- 
figuration p = (0, 0, 0) in the following. The remaining dimensions are only the 
time direction and spinor space, therefore for a given Lg and the fixed spatial 
momentum configuration, A^-^^^ is a matrix with dimension 4{T/a — 1). Figure 
[T] shows the magnitude of ln(||A^^''-'(a;o,yo)||spin)) where the norm is taken for 
the spinor space only. By increasing Lg, the bulk symmetry breaking is reduced. 
Finally at Lg = 24, only boundary breaking effects remain and they are localized 
exponentially near the time boundaries. This is the expected behavior for overlap 
fermions [12] . We conclude that the presence of the boundary term causes chiral 
symmetry breaking, which decays exponentially away from the time boundaries 
for the effective four dimensional operator. 

3 Spectrum of free operator 

In this section, we investigate the free spectrum of the DWF operator to confirm 
universality at the tree level. We set T = L in this section. 
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c-0,U-4 




- c-0,L-16 




c-0,L3-32 




r 









0.05 0.1 0.15 0.2 

a/L 

Figure 2: The lowest eigenvalue of L^D^^pDowF with am^ = 1. Some combinations 
of parameters c and Lg are shown. The continuum value is 7r^/4 = 2.467... For c = 1 
case, since the Ls-dependence is so weak on on this scale, we show only Lg = 4 results 
as representative. 




Figure 3: ||'?/'(a;o, s)||spin with zero spatial momentum and the parameters T/a = Lg = 
40, ^ = and am^ = 1. The norm for the eigenvector is taken in the spinor space. The 
left (right) panel is for c = 1 (c = 0). 



3.1 Spectrum of D^^^pDowF 

To achieve better chiral symmetry, the (physical) eigenmodes of the domain- 
wall operator should be localized near the boundaries of the fifth direction and 
propagate in the space-time directions. This should also be true in the SF setup. 
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since the chiral symmetry is supposed to be maintained in the bulk. To observe 
such phenomena, we numerically compute the lowest eigenmode of the operator 
-^DWF^DWF with the trivial gauge configuration U{x,fi) = 1. In the free case, 
we can perform the Fourier transformation for spatial directions, and project out 
the momentum configuration p = (0, 0, 0). Thus, remaining indexes of the vector 
space are now the spinor, the time xq and the extra dimension s, 

^DWF-^DWF'0(2;o, s) = Xi^ixo, s), (35) 

where the spinor indexes are suppressed. We set input parameters ams = 1 and 
9 = 0. 9 is the parameter which controls the spatial boundary condition for 
fermion fields. (For more details, we refer to [18].) 

We numerically compute the lowest eigenvalue and the corresponding eigen- 
function ||?/'(xo, s)||spin- We examine not only for c = 1 but also for c = to 
investigate the importance of the presence of the boundary operator. The scal- 
ing behavior of the eigenvalue is shown in Figure [2j For c = 1, Lg dependence is 
too small to see on this scale. The lowest eigenvalues converge to their continuum 
values properly, therefore universality is confirmed. Furthermore, the associated 
eigenfunction shows nice localization behavior, namely being localizing for the 
fifth direction and propagating for the time direction, as shown in the left panel 
of Figure El This shows that this mode is a physical one. 

For c = in Figure [2l although all Lg = 4, 16, 32 results tend to converge 
to the continuum limit, large Lg results have a bending phenomenon in small 
a/L region and show no power decay in terms of a/L. This indicates that if 
one takes Lg to infinite before taking a/L = limit, the eigenvalue will likely 
converge to zero. If this is so, the theory with Lg = oo does not belong to a 
correct universality class. Furthermore, the eigenfunction in the right panel in 
Figure [3] is localized on edges in the time-s plane. This is a typical unphysical 
mode. On the other hand, interestingly for small Lg, the scaling behavior is 
rather mild. In the small Lg case, the chiral symmetry breaking of domain-wall 
fermions are rather similar to that of the ordinary Wilson fermions. As in the 
Wilson fermions case, the bulk chiral symmetry breaking for DWFs due to finite 
Lg plays some role in producing the correct continuum limit. This is the reason 
why DWFs with smaller Lg and no boundary term B can produce the continuum 
results. 

The results shown in this subsection show that the boundary term with c 7^ 
plays an important role for the theory to be in the correct universality class. 

3.2 Spectrum of DjDq 

Not all eigenmodes of -Ddwf are physical ones and elimination of unphysical 
mode is not clear. To extract physical modes only, let us study the eigenmodes 
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I 1 1 1 1 I 1 1 1 1 

0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 

a/L a/L 

Figure 4: The a/L dependence of the lowest ten eigenvalues of L'^D'^^Dq in the presence 
of the background gauge field. The left (right) panel is for 6 = {6 = vr/5). The 
parameters are set to am^ = 1 and c = 1. The red points at a/L = are continuum 
values [T8] . 



of -Dq, where unphysical modes are excluded. The operator Dq is defined from 
Sq in eq. (f33]) 

DqSq{x,y) = a'^dx^y. (36) 

We numerically compute the lowest ten eigenvalues of DqDq with the param- 
eter set oms = 1, 9 = 0,7r/5 and c = 1 in the presence of the the background 
gauge field (choice A in Ref. [IH]). The values obtained for L/a = T/a = Q, 12, 24 
are summarized in Table [TJ All tables are given in appendix [B] The scaling be- 
havior of the eigenvalues are shown in Figure HI Although we show two cases 
of Lg, namely = 4 and Ls = 32, it is hard to see the difference on this 
scale. We observe that they converge to the continuum values given in Ref |18j . 
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This behavior persists for a variety of values of c, 0.5 < c < 1.5. This confirms 
universahty at the tree level. 



4 One-loop analysis of SF coupling 

To check further universality at the quantum level and renormalizability, we 
perform the one-loop order calculation of the SF coupling. 



4.1 Definition and results 



We compute the fermion contribution to the SF coupling [18j pi^i{L/a,Ls) (we 
set L = T as usual) at one-loop order for massless domain-wall fermions. The 
one-loop coefficient is given as 



1 d 

- — In det (Z^DWF / Dpy ] 



pi^i{L/a,Ls) 

with a normalization (See [18] for details.) 

k = 12(L/a)^[sin(7) + sin(27)] 



rj=i/=0 



7 



-7r{a/Lf 



(37) 



(38) 



The parameters tj and v parameterize the background gauge field [19] . In the ac- 
tual calculation, we expand the rj derivative and use the fact that the determinant 
is factorized for individual spatial momentum p and color sector b, 



pi^i{L/a,Ls 



1 



Tr 



D 



1 

DWF' 



1 ^ 



drj 
(-C'dwf 



D 



drj 



PV 



p b=l 



b ^-i^p-j_^-^r)WF(p) 



dri 



-iD'py)-\p) 



drj 



(39) 



The trace tr concerns the spinor, the time indices and fifth coordinate only. It 
is maybe worthwhile to note that for our definition of DWF, 



(9Z)dwf dDpY 



drj 



drj 



(40) 



holds since the mass term does not involve the gauge field. 

We compute pi^i on the lattices of size L/o = 4, 6, ...,48 and Ls = 6, 8, 10, 12, 16 
with parameters 0.7 < am^ < 1.3 and = 7r/5. Subsets of the results are 
summarized in Table [2] for am^ = 1, Lg = 6 and L/a = 4, 6, ...,48. Separate 
contributions from DWF and PV are also shown there. 
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Figure 5: L<j-dependence of Aq for 0.7 < am^ < 1.3. The horizontal hnes show the 
values of Aq in the infinity Ls limit, which are obtained by combining the results of 
previous literature [TBI EO] . 



4.2 Coefficients of Symanzik's expansion 

Prom the Symanzik's analysis of the cutoff dependence of Feynman diagrams on 
the lattice, one expects that the one-loop coefficient has an asymptotic expansion 
in terms of a/L 

oo 

Pi,i{L/a, Ls) = J2{a/Lr[An{Ls) + B„(L,) ln(L/a)]. (41) 

n=0 

Note that the coefficients An and Bn (n = 0,1,2,...) depend on Ls- We can 
reliably extract the first few coefficients by making use of the method described 
in Ref. [21] ■ 

For the usual renormalization of the coupling constant, Bq at Ls = oo should 
be 26o,i where 6o,i is the fermion part of the one- loop coefficient of the /3-function 
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for Nf flavors QCD, 



bo = 
bo,o = 

bo,i = 

We confirmed that Bq[Ls) for large Ls (say = 16) converges to 26o,i = 
—0.008443... up to three significant digits for the values of am^ which we in- 
vestigated. When the tree-level 0(a) improvement is realized, we expect that 
5i = holds. We check this to 10"^ for the same parameter region as before. 
This shows that the formula for the boundary coefficient in eq. (j27p works well 
to achieve the tree- level 0(a) improvement to the precision considered here. In 
the following analysis, we set exact values Bq = 26o,i = ~l/(12vr^) and Bi = 0. 

Aq gives information about a ratio of A-parameters. The obtained values of 
Aq{Ls) as a function of Lg are shown in Figure El By combining the previous 
results from Ref. [18tl20|. the values of Aq at infinity Lg can be obtained, and are 
shown in Figure [5] as the horizontal lines. We observe that our results at finite 
Lg properly converge to the known results at infinity Lg. 

To achieve one-loop 0(a) improvement, we need to determine the coefficient 
of the fermion part of the boundary counter-term, c|^'^'' [18j at one-loop order. 
If one imposes an improvement condition |18j . one finds that 

''^ = Ai/2, (45) 

therefore we need the value of Ai. The obtained values of Ai are given in Tabled 
For future reference, we provide an interpolation formula for c^'^^ as a polynomial 
of ams for larger L^, where value of Ai is saturated, 

c|^'^) = 0.00434 + 0.01102(am5 - 1) - 0.00858(am5 - 1)^, (46) 

for 0.7 < ams < 1.3. 

5 Lattice artifacts of the step scaling func- 
tion to one-loop order 

In this section, we investigate lattice artifacts of the step scaling function (SSF) 
[22] a{2,u), which describes the evolution of the running coupling (^{L) = u 
under changes of scale L by a factor 2, 

a{2,u)=f{2L), u=-g\L). (47) 



00,0 + ^^fO0,l, 
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(43) 

(44) 
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Figure 6: We show the relative deviation with the various actions for tree level 0{a) 
improvement, 6^1 (Left), and one-loop 0{a) improvement, S^i^l (Right), as a function 
of a/L and {a/Ly respectively. Upper part is for 6* = 0, and lower is for ^ = vr/S. For 
comparison, those of the Wilson type fermion with 4^'^"* = and the clover fermion 
with cj^'^^ = 0.019141 [18] are included in the plot of Sf'l and 6^1 respectively. 



The lattice version of the step scaling function is denoted by S(2,n, a/L) which 
contains lattice artifacts. Such lattice artifacts are described by the relative 
deviation 

_ S(2,u,a/L) - cr(2,n) 
d{u,a/L) = — . (48) 

(T{Z, U) 

By expanding the relative deviation in terms of the coupling constant u, one 
obtains 

6{u,a/L) = 6i{a/L)u + 0{u^), (49) 

where the one-loop deviation, 6i{s,a/L), may be decomposed into pure gauge 
and fermion part |18] . 

6i{a/L) = 5i,o(a/L) + iVf Ji,i(a/L). (50) 

We are currently only interested in the fermion part. We consider domain- wall 
fermions, thus the fermion part of the one-loop deviation 6i^i{a/L, Ls) contains 
Lg dependence. In terms of the one-loop coefficient of the SF coupling pi^i, the 
one-loop deviation is given by 

5i,i(a/L, Ls) = Pi,i(2L/a, L,) - pi,i(L/a, L,) - 26o,i ln(2). (51) 

Depending on the value of the boundary counter term c{:^'^\ we denote with 6^1 
the tree level 0(a) improved version with cf '^^ = 0, and 5^^l the one- loop 0(a) 
improved one for c^'^^ in eq. (j45p . 
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We show numerical results for the one-loop deviation in Tables H] and [5] and 
the plots in Figure [U where we include those of the Wilson type fermions |18j 
and overlap fermion [13j for comparison. Lg-dependence of DWFs is small. In 
the case of the clover action, c[^'^^ is set to be the proper value to achieve one- 
loop 0(o) improvement, and for Wilson fermions it is set to c[^'^^ = 0. We 
observe that the lattice artifacts for domain-wall fermions are small for tree level 
boundary 0(a) improvement case compared with other fermions, while they are 
large for one-loop boundary 0(a) improvement case. 

6 Conclusion and outlook 

In this paper, we provide a new formulation of domain-wall fermions in the SF 
setup by following the universality argument of Liischer. In contrast to the pre- 
vious formulation by Taniguchi, ours can deal with the boundary 0(a) improve- 
ment properly, and there is no constraint on the number of flavors. To check that 
our formulation works properly, we investigate the spectrum and eigenmodes of 
the free operator, and perform a one-loop analysis of the SF coupling constant. 
Then we confirm universality at tree and the one-loop level and observe that all 
results investigated show the desired behaviors. 

Before starting simulations, the boundary improvement coefficient c should 
be determined to one-loop order. This involves calculations of the SF correlators, 
/A) /p etc. given in appendixO This could be done in a similar way to the case 
of the Wilson fermion. 

As mentioned before, one of the most important properties of the universality 
formulation is that there are no restriction of the number of flavors. By taking 
advantage of this property, we may compute the renormalization factor of 
for A^f = 3 QCD. 
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A Fermion correlators 



In this appendix, I summarize fermion correlators, the boundary fields and so on 
which are often used in the SF setup. 

A.l Boundary fields 

As given in Ref. p^, the lattice version of the fermion boundary fields are defined 

C(x) = U{x-a6,0)P^q{x)\„, (52) 

C(x) = q{x)P+U{x-aO,0)-\,=a, (53) 

C'(x) = U{x,Or'P+q{x)Uo=T-a, (54) 

C'(x) = q{x)P^U{x,0)U,=T-a. (55) 



Here note that we use the physical quark fields defined in eq. (pO|) and (j3T|) . 

A. 2 Propagators 

The propagators for the physical quark fields and the boundary fields are given 
by 

[Qi^My)]F = Sq{x,y), (56) 

[9(^)C(y)]p = S^{x,y)Uiy-a6,0)-^P+\y,=a, (57) 

['7(a;)C'(y)]p = S^{x,y)U{y,0)P^\y,=T-a, (58) 

[C(x)g(y)]p = P^U{x-a6,0)S^{x,y)\„, (59) 

[C'(x)^(y)]p = P+[/(x,0)-i5q(x,y)Uo=T-a, (60) 

[C(x)C'(y)]p = P.U{x-a6,0)S^{x,y)U{y,0)P.Uo=a,y,=T-a, (61) 
[C'(x)C(y)]p = P+Uix,0)-^S^{x,y)Uiy - a6,0)-^P.\^,=T-a,yo=a. (62) 

A. 3 Operators 

We consider the degenerate quark mass case and an extension of the flavor space 
is done in a trivial way. In terms of the physical quark fields, the local operators 
are defined as 

Al.ix) = q{x)-f^j5lr\{x), (63) 
P"(x) = q{xh,Ky{x). (64) 
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The conserved axial vector current is given by 



= ^ sign Is 



s=l 



Ls + l 



2 r'' 



where 



f^{x,s) = tPix + afi,s)Pi^^U{x,^l)-'^T^^l;{x,s) 
-i;{x,s)p[^^U{x,n)-T'''i(j{x + ail,s), 



withPi'^) = (1±7m)/2- 

A. 4 Correlators 

The fermion correlators for the local operators in the SF are given by 



/A(a;o) 
fp{xo) 



EE]^^^^o(^)C(y)75 2T'^C(z)), 

a=l y,z f 

Z^ Z^ N2 

a=l y.z I 



:^)C(y)75ir"C(z)), 



.12 



/i = "76 E E ]^2^(C'(u)75^r'^C'(v)C(y)75^r'^C(z)). 



a=l u,v,y,z " f 

After Wick contraction, they become 

/a(xo) = a6 5;i([C(z)g-(x)]F7075[g(x)C(y)]75), 

/p(xo) = a6^i([C(z)g-(x)]F75k(x)C(y)]75), 

/i = E J([C(^)C'(u)]F75[C'(v)C(y)]F75), 



u,v,y,z 



where the propagators are given in subsection IA.2I and we have used 

2' 



a=l 



Nf-1 



These correlators are the same as those of Wilson fermions except that the 
agators are replaced by those of the physical quark field Sq. 
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For the conserved axial vector current, a correlator is given by 



fA{xo) 



1 



(-4g(x)C(y)75^T'^C(z)). 



(74) 



N?-l 



a=l y,z 



As an example, at tree level, this can be expressed in terms of the propagator 
for domain- wall fermions as. 
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e 


= 












L/a 


= 6 


L/a 


= 12 


L/a 


= 24 






n 


Ls=4: 


Ls = 32 


Ls = 4. 


Ls = 32 


Ls=4: 


Ls = 32 


b 


d 


1 


3.161141 


3.160760 


2.591269 


2.591267 


2.350053 


2.350053 


2 


2 


2 


5.658392 


5.658925 


5.191148 


5.191161 


4.990293 


4.990293 


2 


2 


3 


9.050173 


9.045312 


8.196497 


8.196424 


7.888424 


7.888424 


3 


2 


4 


11.981137 


11.963010 


10.635061 


10.634737 


10.177881 


10.177878 


1 


2 


5 


13.098016 


13.101850 


12.434736 


12.434839 


12.281233 


12.281235 


3 


2 


6 


22.037107 


22.078353 


20.612026 


20.613136 


20.436083 


20.436098 


1 


2 


7 


30.304338 


30.232012 


25.944369 


25.942289 


24.578258 


24.578235 


2 


2 


8 


30.810696 


30.865401 


26.708727 


26.709232 


25.378691 


25.378700 


2 


2 


9 


30.585955 


30.530965 


27.241272 


27.239584 


27.184415 


27.184395 


1 


6 


10 


31.026393 


30.978110 


28.357563 


28.355829 


28.438888 


28.438866 


3 


6 








e = 


= 7r/5 












L/a 


= 6 


L/a 


= 12 


L/a 


= 24 






n 


Ls=4: 


Ls = 32 


Ls=4: 


L, = 32 


Ls=4: 


Ls = 32 


b 


d 


1 


5.924559 


5.922886 


5.232916 


5.232896 


4.952553 


4.952553 


2 


2 


2 


6.428276 


6.423570 


5.566868 


5.566810 


5.214696 


5.214695 


1 


2 


3 


9.221621 


9.223322 


8.721989 


8.722031 


8.548948 


8.548949 


2 


2 


4 


13.912656 


13.926143 


13.392223 


13.392552 


13.267533 


13.267537 


1 


2 


5 


15.970852 


15.951467 


14.513964 


14.513548 


14.162258 


14.162253 


3 


2 


6 


21.508691 


21.520321 


19.983404 


19.983750 


19.838031 


19.838036 


3 


2 


7 


35.296282 


35.214384 


29.231042 


29.228128 


27.743584 


27.743548 


2 


2 


8 


35.861910 


35.932144 


30.057783 


30.058382 


28.612335 


28.612346 


2 


2 


9 


30.468452 


30.406862 


27.727188 


27.725168 


27.806412 


27.806386 


1 


6 


10 


31.520474 


31.475867 


28.021244 


28.019681 


27.896264 


27.896245 


3 


6 



Table 1: The lowest ten eigenvalues of the Hermitian operator L^Dj^D^ for Lg = 4, 32. 
Upper (Lower) panel is for 6* = (6* = 7r/5). b represents the color sector, and d is for 
degeneracy for one flavor. 
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L/a 


pi,i(L/a,6) 


DWF contribution 


PV contribution 


4 


-0 


0090558230 


-0 


0463742376 


-0 


0373184146 


6 


-0 


0128831139 


-0 
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-0 


0485656342 
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-0 
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-0 


0835634805 
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-0 
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-0 
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-0 


0750539311 


44 


-0 


0314533516 


-0 
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-0 


0752530061 


46 


-0 


0318432669 


-0 


1072780356 


-0 


0754347687 


48 


-0 


0322159188 


-0 


1078173020 


-0 


0756013832 



Table 2: The one-loop coefficient of the SF coupling pi i(L/a, Lg) with = 6, am^ = 1 
and 9 = tt/5. In eq. fl39|) . there are two sources of contributions: DWF and PV, which 
are shown separately in the table. 



Ls \ am^ 


0.7 


0.8 


0.9 


1.0 


1.1 


1.2 


1.3 


6 






0.0102(8) 


0.0125(9) 


0.0145(7) 






8 




0.0047(4) 


0.0074(9) 


0.0097(9) 


0.0119(9) 


0.0135(4) 




10 


0.0004(2) 


0.0040(9) 


0.0066(9) 


0.0090(9) 


0.0111(9) 


0.0129(9) 


0.0135(2) 


12 


0.0007(9) 


0.0037(9) 


0.0064(9) 


0.0088(9) 


0.0108(9) 


0.0126(9) 


0.0139(8) 


16 


0.0006(10) 


0.0036(9) 


0.0063(9) 


0.0087(9) 


0.0107(9) 


0.0125(9) 


0.0137(10) 



Table 3: The value of Ai for Lg = 6, 8, 10, 12, 16 and 0.7 < am^ < 1.3. 
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L/a \ Ls 
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-0.000545 
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-0 
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-0.000284 
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24 


-0.000322 


-0.000273 


-0 
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-0.000256 


-0.000254 



Table 4: The relative deviation 5{^{ with am^ — 1 and 6* = 7r/5 for tree level boundary 
0(a) improvement. 



L/a \ Ls 


6 


8 


10 


12 
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0.000955 


0.000701 


0.000597 


0.000558 


0.000539 


6 
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10 
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12 
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-0.000194 
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14 
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-0.000162 


-0.000166 
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16 


-0.000123 


-0.000134 


-0.000138 


-0.000139 


-0.000140 


18 


-0.000101 


-0.000112 


-0.000115 


-0.000117 


-0.000117 


20 


-0.000085 


-0.000094 


-0.000098 


-0.000099 


-0.000099 


22 


-0.000072 


-0.000080 


-0.000083 


-0.000084 


-0.000085 


24 


-0.000062 


-0.000069 


-0.000072 


-0.000073 


-0.000074 



Table 5: The relative deviation 6^1 with am^ — 1 and 6 = n/b for one- loop level 
boundary 0(a) improvement. 
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